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ABSTRACT 


The  report  investigates  scattering  from  a  rough  surface  when  trans¬ 
mitter  and  receiver  are  at  a  finite  distance  from  each  other,  so  that 
the  sphericity  of  the  incident  and  scattered  waves  has  to  be  taken  into 
account.  The  calculation  is  complicated  by  the  fact  that,  unlike  the 
case  of  plane  waves,  the  equiprobability  planes  of  the  surface  do  not 
coincide  with  the  equiphase  surfaces,  and  therefore  a  number  of  simpli¬ 
fying  assumptions  are  made,  in  particular,  geometric  optics  are  used. 

The  phase  distribution  of  the  scattered  field  by  approximate  Integra- 

A 

tions.  From  this,  the  other  quantities  of  interest  can  be  found.  The 
results  are  compared  with  the  measurements  performed  by  Beard  using 
microwaves  and  a  random  water  surface  with  known  statistical  properties. 
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1.  INTRODUCTION 


When  a  rough  surface  generated  by  a  stationary  random  process  is 
illuminated  by  a  plane  wave  and  the  scattered  field  is  observed  far  in 
the  Fraunhofer  zone  (i.e.,  the  scattered  waves  may  be  considered  plane) 
the  reflection  points  on  the  surface  resulting  in  the  same  phase  $  at 
the  receiving  point  are  on  a  plane  z  *  const,  where  z  is  measured  from 
the  mean  plane  of  the  surface.  These  equiphase  planes  are  thus  identi¬ 
cal  with  the  equiprobability  planes  Z  *  z,  where  Z  is  the  deviation  of 
the  surface  from  its  mean  plane,  a  random  variable  (here  assumed  normal) 

The  ultimate  root  of  the  difficulties  encountered  in  an  investiga¬ 
tion  of  non-plane  wave  illumination  and  scattering  is  the  fact  that  the 
equiphase  surfaces  are  no  longer  identical  with  the  equiprobability  sur¬ 
faces.  The  equiprobability  surfaces  are  given  by  p(z)  *  const,  where 
p(z)  is  the  probability  density  of  Z;  as  before,  the  equiprobability 
surfaces  are  therefore  given  by  the  planes  z  *  const  (Figure  1).  How¬ 
ever,  the  locus  of  reflection  points  yielding  the  same  phase  at  the  re¬ 
ceiver,  when  both  are  at  a  height  h  above  <Z>  and  separated  by  a  dis¬ 
tance  D,  is  a  rotational  ellipsoid  with  foci  at  the  transmitting  and 
receiving  points,  T^  and  R^,  respectively. 

The  resulting  phase  distribution  at  the  receiver  will  under  these 
conditions  be  asymmetrical,  which  leads  to  further  consequences  rarely 
met  in  plane  wave  propagation.  To  see  qualitatively  why  the  phase  dis¬ 
tribution  is  asymmetrical,  consider  the  case  when  the  standard  deviation 
ow  of  the  surface  is  very  small  so  that  the  surface  deviates  only  very 
slightly  from  a  plane  (Figure  2).  For  plane  wave  illumination  and  scat¬ 
tering,  this  would  obviously  result  in  a  received  phase  that  only  very 
slightly  deviates  from  its  mean.  However,  for  spherical  waves,  the 
phase  is  given  by  4>  *  k(rj  +  r2).  Let  the  point  x  *  0,  y  ■  0  (when 
Z  ■  0)  give  rise  to  a  phase  $0.  Smaller  phases  are  possible,  though  not 
very  probable;  they  will  occur  only  when  Z(0,0)  goes  positive.  On  the 
other  hand,  phases  larger  than  4>0  will  occur  very  often,  since  practi¬ 
cally  any  reflecting  point  with  coordinates  other  than  x  *  0,  y  «  0  will 
give  rise  to  a  phase  iarger  than  #0.  The  phases  near  the  value  $  will 


2 


I 

. . . 


occur  for  reflecting  points  in  the  first  Fresnel  zone,  and  since  this  is 
the  largest,  the  probability  density  of  $  can  be  expected  to  have  a  max¬ 
imum  near  this  value.  For  the  reasons  just  given,  this  probability  den¬ 
sity  should  then  fall  off  quickly  from  this  maximum  towards  smaller  <t>, 
and  slowly  towards  larger  $  as  shown  in  Figure  3.  Thus,  p($)  will  be 
asymmetrical. 

The  main  effort  of  this  report  is  to  derive  this  phase  distribu¬ 
tion;  once  this  is  established,  the  other  quantities  of  Interest,  in 
partlculer,  the  distribution  of  the  scattered  field,  and  hence  the  mean 
scattered  power,  the  behavior  of  the  phase-quadrature  components,  etc., 
can  be  derived  from  this  phase  distribution  by  previously  established 
methods  (Beckmann,  1962;  Beckmann  and  Splzzichino,  1963). 

The  calculation  is  Involved  and  it  was  not  found  possible  to  arrive 
at  practically  useful  results  unless  a  number  of  simplifying  assumptions 
were  made.  These  are  the  following: 

(A)  Geometric  optics  are  valid;  only  waves  reflected  from  a  favorable, 
or  nearly  favorable,  slope  (making  the  local  angle  of  incidence 
equal  to  the  local  angle  of  reflection)  will  reach  the  receiver. 

(B)  The  equlphase  surfaces  are  elliptical  cylinders  rather  than  rota¬ 
tional  ellipsoids;  near  the  plane  y  ■  0  (where  most  reflections 
must  be  expected)  the  two  are  identical  and  thus  the  problem  is 
reduced  from  a  3-dimensional  one  to  a  2-dimensional  one. 

(C)  The  transmitting  and  receiving  antennas  are  at  the  same  height  h 
above  the  mean  rough  surface. 

(D)  The  rough  surface  is  generated  by  a  stationary  normal  process  in 
the  x  direction. 

Assumptions  (C)  and  (0)  may  be  withdrawn  at  the  cost  of  mathemati¬ 
cal  simplicity,  and  the  same  is  essentially  true  of  Assumption  (B), 
though  this  will  lead  to  integrals  of  almost  prohibitive  complexity;  but 
Assumption  (A)  is  essential  to  the  method  to  be  used. 


Fig.  1.  Equiprobability  and  equiphase  surfaces. 


Fig.  2.  Qualitative  derivation  of  the  phase  distribution 
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In  addition  to  the  above  simplifying  assumptions,  approximate  pro¬ 
cedures  are  used  to  evaluate  some  otherwise  intractable  integrals.  In 
spite  of  this  coarse  treatment,  the  final  results  appear  to  be  at  least 
in  qualitative  agreement  with  experimental  results,  so  that  the  method 
should  be  useful  for  providing  physical  insight,  if  not  precise  numeri¬ 
cal  results. 


2.  PHASE  DISTRIBUTION  OF  THE  SCATTERED  FIELD 


Consider  the  distance  (or  "geometric  nhaae")  *  from  T^  to  R^  via  a 
reflection  ooint  P  (Figure  A).*  For  a  giver  *,  the  reflection  point  must 
lie  on  an  ellipse  with  semimajor  axis 

a  -  is*  (1) 


and  semiminor  axis 


(2) 


where  D  ■  2d  is  the  distance  between  transmitter  and  receiver, 
equation  of  the  ellipse  is 


The 


(3) 


In  the  following  only  the  lower  left  quadrant  of  the  ellipse  is 
considered:  the  upper  quadrants  are  irrelevant  because  the  surface  Z(x) 
has  such  a  small  standard  deviation  that  it  is  practically  impossible 
for  Z(x)  i  h  to  be  satisfied,  and  the  lower  right  quadrant  only  rein¬ 
forces  all  effects  by  symmetry. 


Consider  now  a  fixed  propagation  distance  or  geometric  phase  * 
assuming  that  the  reflection  point  was  located  near  a  height  Z  *  z. 
This  uniquely  determines  the  x  coordinate  of  the  reflection  point  and 
is  determined  from  (3)  for  Z  ■  z: 


x 

z 


(A) 


Upper  and  lower  case  *  and  *,  respectively  the  geometric  and  elec¬ 
trical  nhases,  are  sufficiently  similar  in  form  that  care  should  be 
exercised  to  prevent  confusion. 


The  slope  of  the  ellipse  at  any  point  is  from  (3) 


Z’(x) 


bx 

a/a2  -  x2 


(5) 


and  hence  the  slope  at  the  reflection  point  (if  reflection  towards  the 
receiver  is  to  occur)  must  be,  as  obtained  by  substituting  (4)  for  x, 


Z’(xz) 


a 

z 


(6) 


The  probability  of  this  being  the  case  is  found  from  the  distribu¬ 
tion  of  Z',  which  is  normal  and  stationary,  since  Z  is  normal  and  sta¬ 
tionary,  with  mean  square  value 


o'2  -  <Z'2>  -  -  B"(0) 


where  B(t)  is  the  correlation  function  of  Z(x). 


Thus,  Lhe  density  of  Z'  is 


p(z') 


T  /  x2z'2\ 

71 zf  exp  "  17^ 

W  \  w  / 


where 


w 


(7) 


(8) 


(9) 


This  quantity,  which  has  the  dimension  of  a  length,  is  in  effect  a 
decorrelation  distance;  if  the  correlation  function  is  Gaussian,  t  is 
equal  to  the  correlation  distance  defined  as  the  distance  for  which  the 
autocorrelation  coefficient  decreases  to  e_1 ;  for  other  correlation 
functions,  x  differs  from  the  correlation  distance  only  by  a  constant 
factor.  Note  that  o'  is  the  rms  slope  of  the  surface. 

From  (6),  (7)  and  (8)  we  obtain  the  conditional  probability  density 
of  4>  given  that  the  reflection  takes  place  at  height  z: 
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, .  I  v  t  r  b2 (b2  —  z2)t2  1 

(4|2)  ■  Tm  eltp["  2aV0  *  J 

w  w 


(10) 


where  3>  is  contained  implicitly  in  the  quantities  a  and  through  (1) 
and  (2). 

The  probability  distribution  p(4>)  is  now  obtained  from  the  Theorem 
of  Total  Probability  by  integrating  over  all  z: 


00 

p($)  =  /  p($| z)p(z)dz 

—  00 


(11) 


where  the  density  of  z  is 


p(z) 


2o 


w 


(12) 


Substituting  (10)  and  (12)  in  (11)  we  obtain 


p(4) 


00 

w  -00 


b2T2(b2  -  z2)  (z  -  h) 2 

2a^z2o  ^  2o  2 

w  w 


dz 


(13) 


To  evaluate  this  integral  exactly  is  a  fairly  hopeless  task,  and 
therefore  the  Laplace  method  of  integration  is  used.  We  denote  the  ex¬ 
ponent  by  E(z)  and  find  the  point  zQ  for  which  it  is  maximum  from  the 
condition 


E ' (z)  -  0 


(14) 


Expanding  the  exponent  in  the  neighborhood  of  this  point  (from  which  the 
integral  obtains  its  greatest  contribution),  we  obtain  on  retaining 
terms  up  to  the  quadratic  term. 


p(<J>) 


TeE(Zo)  -  [  E"(z  > 

I  — J - 

W  -O0  L 


(z  -  zo) 


2 


dz 


(15) 
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To  solve  the  cubic  equation  resulting  from  (14),  we  obtain  by  New¬ 
ton's  method 


z  es  h 
o 


(16) 


so  that  after  some  algebra  one  obtains 

T  b2r2(b2  -  h2) 

T  expL'  “Ya'fe- 2 


p(4>) 


s  (<J>) 


where 


’<*>  ■  +  ^r) 


(17) 


(18) 


Substituting  from  (1)  and  (2) ,  this  yields 

p($) 


77(40  exp  ‘ 

w  L 


(<fr2  -  p2)($2  -  p2  -  4h2)t2 
8$2h2a_  2 


w 


(19) 


Note  that  this  expression  has  been  obtained  from  (13)  by  consider¬ 
ing  the  contribution  from  the  region  near  z  *  h,  but  that  the  result 
will  be  accurate  only  if  the  factor  multiplying  (z  -  z0)2  is  large (guar¬ 
anteeing  that  almost  all  of  the  contribution  to  the  integral  is  concen¬ 
trated  into  the  considered  region),  i.e.. 


b2  t  2 
2a2o  2 


>> 


1 


(20) 


In  practical  applications  ow  <<  t,  so  that  (20)  will  be  well  satis¬ 
fied  for  probable  phases  4>,  for  which  a  is  of  the  same  order  as  lb.  How¬ 
ever,  for  the  highly  improbable  phases  $  -*■  D  (corresponding  to  the  sur¬ 
face  rising  to  a  height  where  it  will  interrupt  the  direct  ray  from 
transmitter  to  receiver),  the  ellipse  will  approach  the  segment  TxRx 
b  -*■  0,  and  (20)  is  violated.  Thus  (19)  is  valid  only  for  probable 
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phases  or  larger.  To  restore  the  validity  of  (19)  we  therefore  only 
need  to  impose  a  lower  limit  on  <J>.  This  is  found  as  follows.  We  con¬ 
sider  the  normal  distribution  of  the  surface  truncated  at  h  -  3ow  (in¬ 
curring  an  error  of  less  than  0.5  percent),  which  corresponds  to  a  mini¬ 
mum  phase  (attained  only  at  the  point  x  ■  0) 

i>  =  /D2+  4(h  -  3o  )2  (21) 

m  w 


or,  since  by  assumption  °w  <<  h* 

4>  -  /d2  +  4h2 

m 


(22) 


The  phase  distribution  may  then  be  written  in  the  form 


p($) 


Ct 


o  s  (4>) 
w 


exp  - 


($2  _  d2)(<j>2  -  <J>  2jt2 
_  m 

85>2h2c 


w 


(23) 


where  the  normalization  constant  C  is  found  from  the  condition  that  (23) 
must  integrate  to  unity.  To  perform  this  integration,  the  LaDlace 
method  of  integration  is  again  used.  The  exponent  is  now  maximum  for 

$  =  /$~D  (24) 

o  m 


But,  since  from  (22) 

♦.  -  D(l  +  w)  (25> 

and  in  practical  applications  h/D  <<  1,  we  are  justified  in  writing  (24) 
as 


$  _  $ 
o  ~  m 


(26) 


Substituting  for  s($)  from  (18)  and  noting  that 

(<i>  2  -  D2)2 

b4($  ) - — T7 - h4  (27) 

ra  lb 
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we  find 


s(4>  ) 
m 


-  V2-(i  +  dTtW)W^ 


(28) 


so  that  the  phase  distribution  is  from  (23) 


j-J5 [ 

hO  "  TT 


P($)  ’  ^ 
w 


t2($2  -  $  )2(4>2  -  D2) 

- 5W7! - 1,2/21)2 


W 


(29) 


or»  if  it  is  approximated  by  its  expansion  near  its  maximum  as  above, 

t  V)  f  t2(4>  -  )2  "1 

p<4)  •  STY?  “>  -  2h20  2 -  (30) 


Transforming  from  the  geometric  phase 
through 


2  ir 

$  *  k$  *  — 


4>  to  the  electrical  phase 


(31) 


we  finally  find 


p(4>) 


t  j;  r  (*  ~  y2,2l 

kho  Vx  exp  2k2h2o  ^ 
w  L  w  J 


($<({>£  00) 
m 


(32) 


where 


<j>  ■  k$  ■  k^D2  +  4h^ 

m  m 


(33) 


This  is  the  phase  distribution  to  be  used  for  the  calculation  of 
the  scattered  field. 


3.  ROTATION  OF  THE  VARIANCE  ELLIPSE 

The  phase  of  the  received  field  is  a  random  variable,  and  the  re¬ 
ceived  field  is  therefore  conveniently  decomposed  into  two  components; 


...  *«A%xnf 


} 
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one  in  phase  with  the  field  reflected  by  a  smooth  plane,  and  the  other 
perpendicular  to  it.  These  two  components  will  be  referred  to  as  the 
real  and  imaginary  component,  respectively,  of  the  scattered  field, 
since  it  is  convenient  to  plot  the  received  field  in  the  complex  plane. 
Since  both  components  represent  the  sum  of  a  large  number  of  random 
quantities,  namely  the  waves  scattered  from  each  of  the  reflection 
points  on  the  surface,  both  components  will  be  normal  and  the  received 
field  is  therefore  described  by  a  two-dimensional  normal  distribution 
(real  component  X,  imaginary  component  Y) .  The  equiprobability  curves 
for  the  two-dimensional  normal  distribution,  as  is  easily  shown  by  ana¬ 
lytical  geometry,  are  ellipses  centered  at  X  ■  <X>,  Y  ■  <Y>;  their  axes 
have  identical  directions,  one  of  which  makes  an  angle  ?0  with  the  real 
axis  X  (the  other  direction  differs  by  90°).  The  angle  ?0  is  found  from 
the  relation 


tan  2? 

o 


2  coy  (X,Y) 
0 „2  -  0  2 


(34) 


where  ax2,  aY2  are  the  variances  of  the  X  and  Y  components,  respectively. 
Thus,  the  ellipses  of  the  family  of  equiprobability  curves  differ  only 
by  their  major  and  minor  axes;  the  centers  and  directions  are  identical. 
One  oi  these  is  the  "variance  ellipse";  this  is  the  one  whose  major  and 
minor  axes  oel,  ae2  (not  necessarily  respectively),  are  related  to  the 
variances  of  X  and  Y  through 


el 


+  o 


e2 


+  o. 


(35) 


It  can  be  shown  that 

ov2  ■  o  2  cos2?  +o  2  sin2? 

X  el  o  e2  o 

a  2  ri  o  2  sin2?  +o  2  cos2? 

Y  el  o  e2  o 


(36) 


and  hence 
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o  2  cos2£  -  a  2  sin2S 
2  X  o _ Y  o 

ei  cos2£ 

o 

o  2  cos2C  -  o  2  sin2£ 

2  Y _ O _ X _ O 

°e2  cos2' 


(37) 


To  relate  these  quantities  to  the  phase  distribution  p(4>),  we  note 
that  in  our  geometric-optics  model  the  received  field  is  produced  by  the 
interference  of  the  waves  reflected  by  n  reflection  points  on  the  sur¬ 
face,  so  thac  (apart  from  an  unknown  constant  of  proportionality  which 
can  be  absorbed  into  the  equally  unknown  number  n) 

n  n 

X  «*  l  cos  $  Y  “  l  sin<*>  (38) 

J-l  J  J-l  J 


Denoting  the  means  of  X  and  Y  by 

a  *  n Jp(4>)cos<f>d4> 

6  *  n  /p  ( d>) sin 4><i4> 

we  find  the  variances  in  the  form 

)cos2<td  6 


j^l  +  |p(<j>)cos2(td 6  -  ^-J 
y|l  -  fp(<t>)cos2<pd<t>  - 


(39) 

(40) 


(41) 

(42) 


The  covariance  of  X  and  Y  is  after  a  little  manipulation  using  the 
fact  that  the  4> j  are  independent  found  to  be 

cov(X,Y)  =  n/p(4>) cos 4>sind>d<b  -  (43) 

so  that  on  substituting  the  last  three  results  in  (34)  we  obtain 
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tan2c  .  fpMsift.ZMl  -■  -ffiZ-n.2T 

0  /p ( 4>>cos24»d^>  -  -  ~7~^— 


(44) 


This  relation  reveals  a  marked  difference  between  the  scattering  of 
plane  and  spherical  waves  from  a  normal  (or  any  other  symmetrically  dis¬ 
tributed)  surface:  For  plane  waves,  the  phase  distribution  of  a  symmet¬ 
rically  distributed  surface  is  also  symmetrical,  so  that  the  denominator 
of  (44)  vanishes  (or,  if  other  coordinates  are  introduced,  is  constant); 
for  spherical  waves  emanating  from  the  transmitter  at  a  finite  distance 
from  the  surface,  however,  we  found  a  symmetrically  (normally)  distri¬ 
buted  surface  to  result  in  an  asymmetrical  phase  distribution  (32),  so 
that  as  the  variance  of  the  surface,  and  hence  the  variance  of  $  in¬ 
creases,  £0  will  by  (44)  change,  i.e,,  the  variance  ellipse  will  rotate. 

Relations  (39)  through  (44)  establish  the  connection  between  the 
scattered  field  and  the  phase  distribution.  In  the  next  section,  these 
quantities  will  be  evaluated  when  (32)  is  substituted  for  p(4>). 

4.  CALCULATION  OF  THE  PHASE  QUADRATURE  COMPONENTS  OF  THE  SCATTERED  FIELD 

To  find  the  mean  of  the  real  component  of  the  scattered  field,  we 
substitute  (32)  in  (39);  on  setting  $  “  ■  t,  this  becomes 

“ "  /  expr~  2k^^-z]cos(t + *„>dt  <45) 

w  o  L  w  J 


On  expanding  the  cosine,  the  two  resulting  integrals  can  be  evalu¬ 
ated  by  using  the  integrals 


j  e'At2cosBtdt  -  e'B2/4A  (46) 


and 


00  rs 

e  -A 

J  e  sinBtdt 
0 


A  e-»2/4A  p  (1.1.82/4*) 

2a  ri\2’2’  ' 


(47) 
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where  jFj  is  the  confluent  hypergeoraetric  function.  After  some  tedious, 
though  straightforward,  calculations,  the  result  can  be  expressed  in 
terms  of  the  parameter 


k2h2o  2 

w 


(48) 


and  the  auxiliary  variables 


-K 

e  , 


n  ■ 


-K 

e 


,  (LI. 

1  \ 2*2* 


K) 


(49) 


Then  (45)  yields 


a  -  n(£cos4>  -  rising  )  (50) 

m  m 


Similarly,  for  the  mean  of  the  imaginary  component  one  finds 

6  ■  n(£sin<J>  +  ncos<J>  )  (51) 

m  m 


To  evaluate  the  integrals  occurring  in  (44),  we  introduce  the  aux¬ 
iliary  quantities 


-4k 

u  *  e 


v 


-4k 

e 


F  /“•—•4k) 
l' l\2,2*  / 


(52) 


One  then  finds 


/  p ( d>> cos2 4>d4>  *  ycos2<J>  -  vsin2<f> 

4  m  in 


00 

/  p ( 4>) sln2 *  psin24>  +  vcos2i{> 

1  mm 


(53) 


(54) 


so  that  on  substituting  these  results  in  (44)  one  obtains 
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tan2£ 


tan2<J>  + 
m 


v  -  21 


°  1  -  "  ~cl2in  t.n2* 

y  -  £2  +  n  m 


This  can  be  written  in  the  form 


£  *  <t  +  %arctan  — V  ; 

o  m  y  -  £2  +  n 


or  discarding  the  auxiliary  variables  by  (49)  and  (52), 


1Fi  -  iri(ffK) 


♦  +  *,ctm  2^  ‘ 


Since  by  (48)  k  is  proportional  to  ow2,  this  establishes  the  rota¬ 
tion  of  the  variance  ellipse  with  the  variance  ow2  of  the  rough  surface. 

For  practical  applications,  the  approximate  expressions  of  iF*  for 
small  and  large  k  are  listed  here  for  reference  purposes: 


1  +  K 


F  /I.!  \  _  el 

1*  1\2’2  /  2k 


(k  «  1) 


(•c  »  1) 


For  the  variances  of  the  real  and  imaginary  components  one  obtains 
from  (41),  (42),  and  (53),  (54) 

o  2  -  r-  1  +  ycos2$  -  vsin2<J>  -  2(£cos$  -  nsin$  )2  (60) 

AzL  m  m  m  mj 

o  2  -  1  -  ycos2ifi  +  vsin2$  -  2(£sin$  +  ncos<j>  )2  (61) 

Y  *  L  m  m  m  mj 
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The  total  variance  of  the  scattered  field,  which  is  proportional  to 
the  incoherently  scattered  power  (at  the  receiving  point)  is  therefore 

o2*o2+o2*  n(l  -  £2  -  n2)  (62) 

x  y 


As  was  to  be  expected,  this  is  independent  of  The  above  rela¬ 

tion  also  satisfies  the  physical  requirements 


lim  o2  ■  lim  o2  *  0  (63) 

0  -+0  K  +  O 

w 


lim  o2  =  lim  o2  ■  n  (64) 

a  -**>  <-*» 

w 


To  find  the  ratio  of  the  major  to  the  minor  axis  of  the  variance 
ellipse 


K  - 


/  V*1 

°max  __  /°el  \ 
°min  \°e2/ 


(65) 


where  the  exponent  +1  or  -1  is  chosen  so  as  to  make  K2  >1,  we  take  the 
ratio  of  the  two  expressions  in  (37),  obtaining 


whe  re 


k2 


(67) 


or  substituting  from  (60)  and  (61), 

1  -  £z  -  n2  +  (o2  -  £2  +  w)cos2^ra  +  (2£n  -  v)«in2$n 

k  "  l'  -  -  nz  -  (nz  -  +  u)cos2$  -  (2£o  -  v)sin2i 

m  m 


(68) 
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so  that  k2,  and  hence  K2,  depend  on  <J>m,  i.e.,  on  the  geometric  configu¬ 
ration  determined  (with  accuracy  of  a  fraction  of  a  wavelength)  by  D  and 
h  in  (33).  For  ow  -*  °°,  i.e.,  <  -+  (67)  yields  k  -*•  1,  hence 

lim  K2  -  1  (69) 

O  -x» 
w 


The  limit  a  ->0,  i.e.,  k  -*■  0,  is  more  difficult  to  establish, 
w 

since 


lim  £  «  lim  u  ■  1 
k+0  k+O 

lim  r)  *  lim  v  ■  0 

<->■0  k-K) 


(70) 


so  that  for  <  -*■  0,  (68)  becomes  an  indeterminate  expression  of  tvpe  0/0. 
Expanding  £,  n,  y  and  v  in  terms  of  k,  we  find 


k2 


k(2  -  ±)  -  k(2  -  7)cos2»m  -  sin2»m  +  ... 

<(2  -  i)  +  <( 2  -  ~) cos24>m  +  sin2^  +  ... 


(71) 


so  that 


lim  k2 
tc-*-0 


1  -  cos2(J) 
1  +  cos2<(! 


m 


m 


tan2<f> 

m 


(72) 


After  a  similar  procedure  one  obtains  from  (56) 


lim  £  -  *  ±  (m  -  0,1)  (73) 

o  m  L 


Substituting  (72)  and  (73)  in  (66),  we  have 

lim  K2  =  °°  (74) 

a  -+0 
w 
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It  will  later  be  found  that  (69)  agrees  with  experiment,  but  (74) 
does  not.  The  reason  is  that  (69)  can  be  derived  directly  from  (41), 
(42)  and  (67)  for  any  distribution  p(<£);  on  the  other  hand,  the  above 
derivation  shows  that  (74)  depends  very  critically  on  the  distribution 
p((f>),  for  a  minute  change  in  the  latter  will  lead  to  a  completely  dif¬ 
ferent  limit  (74).  For  example,  if  4>m  =  mn  (m  =  1,  2,  ...),  then  the 
indeterminate  expression  0/0  in  (68)  cannot  be  resolved  until  the  quad¬ 
ratic  terms  (i.e.,  terms  in  ow4)  are  taken  into  account.  Obviously,  the 
approximations  used  in  deriving  p($)  cannot  achieve  the  required  accu¬ 
racy  of  up  to  and  including  fourth-order  terms  in 

This  point  is  noted  here  because  (74)  i  the  <nc  result  that  does 
not  agree  qualitatively  with  experimental  ob.  ■  .  ration 

5.  MODIFICATION  FOR  SURFACES  GENERATED  DV  A  NARROW-!  .,.D  PROCESS 

The  underlying  idea  of  a  saddle-pjint  integration  in  physical  op¬ 
tics  to  calculate  the  scattered  field  is  the  summation  of  the  contribu¬ 
tions  to  the  scattered  wave  from  the  neighborhood  of  a  point  with  favor¬ 
able  slope  on  the  surface.  (Saddle-point  methods  in  physical  optics 
thus  represent  a  transition  to  geometric  optics.)  As  the  distance  from 
that  point  increases,  the  slope  will  change  to  a  less  favorable  one, 
which  means  that  the  scattered  wavelets  will  increasingly  be  out  of 
phase  with  the  one  scattered  from  the  favorable  slope,  and  therefore 
both  the  scattered  wave  and  the  integral  will  cease  to  obtain  contribu¬ 
tions  as  the  immediate  neighborhood  of  a  favorable  slope  is  abandoned. 
This  is  schematically  shown  in  Figure  6  for  a  wide-band  random  process. 

However,  for  a  narrow-band  process  the  situation  is  somewhat  dif¬ 
ferent.  A  (mean-square  continuous)  narrow-band  process  consists  of  ran¬ 
dom  oscillations  modulated  by  a  slowly  varying  random  envelope  (Beck¬ 
mann,  1967,  pp.  249-254).  As  shown  in  Figure  7,  several  successive 
oscillations  will  produce  a  succession  of  points  of  the  type  in  Figure  6 
reflecting  favorable  phases. 


2n 


In  calculating  the  field  by  physical  optics,  the  entire  correlation 
function  (which  itself  also  consists  of  oscillations  modulated  by  a 
slowly  varying  envelope)  is  required,  but  this  problem  will  take  care  of 
itself,  since  there  will  be  a  corresponding  succession  of  saddle  points 
whose  contributions  will  add  up  to  yield  the  total  scattered  field. 
However,  since  the  above  analysis  utilizes  not  the  entire  correlation 
function  (as  it  would  in  physical  optics),  but  only  the  RMS  slope  in 
(9),  the  slope  needed  is  not  the  slope  of  the  random  process  itself,  but 
that  of  its  envelope.  The  reason  is  that  in  geometric  optics  one  looks 
for  the  slopes  that  will  concentrate  waves  with  equal  phases;  it  is  ob¬ 
vious  from  Figure  7  that  this  will  be  the  slope  of  the  envelope  of  a 
narrow-band  process,  and  not  the  slope  of  the  narrow-band  process  it¬ 
self  (which  is  much  larger).  It  can  be  shown  [Beckmann,  1967, 
pp.  250-253,  (7-3.8),  (7.3-26),  (7.3-30)]  that  the  correlation  function 

of  this  envelope  is  equal  to  the  envelope  of  the  correlation  function 

B(t)  -  a(T)co8tDoT  (75) 

of  a  narrow-band  process  whose  spectrum  is  symmetrical  about  some  fre¬ 
quency  oj0;  its  rms  slope  is  therefore 


o 


1 

env 


/-  a"(0) 


o 

w 


T 

env 


so  that 


T 

env 


(76) 


(77) 


Since  from  (75) 


a"(0)  -  B"(0)  +  a  2  u  2 

wo 


(78) 


i 

i 


} 


and  B"(0)  <  0,  we  have 
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Fit;.  6.  Single  saddle  point  used  for  a  wide-hand  process. 


Fig.  7.  Several  saddle  points  contribute  to  the  integra¬ 
tion  for  a  narrow-band  process;  this  is  equivalent 
to  using  a  single  saddle  point  of  the  slowly  varying 
envelope  of  the  process  and  its  correlation  function. 
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o 

_ w _ 

env  =  B"  (0)  |  -  a  2  0)  2 
1  'wo 


where 


w 


is  the  autocorrelation  coefficient  of  the  random  surface. 


(79) 


(80) 


6.  COMPARISON  WITH  EXPERIMENTAL  DATA 

There  have  apparently  been  only  two  detailed  experiments  performed 
to  investigate  scattering  from  a  rough  surface  at  a  finite  distance  from 
transmitter  and  receiver,  one  by  Beard  (1965,  1967),  and  the  other  by 
Clarke  and  Hendry  (1964).  Of  these,  however,  only  Beard  has  measured 
the  parameters  determined  directly  by  the  phase,  namely  £0,  oel,  oe2; 
Clarke  and  Hendry  have  only  measured  the  equivalent  of  ox,  oy,  which  is 
ambiguous  for  determining  the  variance  ellipse,  as  can  be  seen  from  Fig¬ 
ure  5,  for  C0  change  with  ow  as  shown  by  (44),  (57),  and  as  also 

follows  from  Beard's  experiments. 

Beard  (1965,  1967)  illuminated  a  rough  water  surface  by  microwaves 
with  wavelength  \  •  2.83  cm.  The  surface  of  water  in  a  tank  was  per¬ 
turbed  by  injecting  air  or  water  into  the  water  tank  (in  a  region  beyond 
the  illuminated  area).  By  varying  the  injection  rate,  the  standard  de¬ 
viation  ow  of  the  water  surface  could  be  varied.  The  transmitter  and 
receiver  antennas  were  separated  by  a  distance  D  ■  4.115  m  (162  inches); 
both  were  situated  at  a  height  h  *  0.869  m  above  the  mean  water  level, 
yielding  a  grazing  angle  ip  *  22.9°  at  the  origin  Z(0,0)  in  Figure  2. 

Thu  distribution  of  the  deviation  of  the  water  surface  Z  was  mea¬ 
sured  by  Beard  (1965)  and  found  normal.  The  correlation  function  of  the 
surface  was  measured  by  Beard  in  1966;  these  data,  partly  not  yet  pub¬ 
lished,  were  at  the  author's  disposal  to  enable  the  above  theory  to  be 
compared  with  experiment.  These  measurements  show  that  the  random 
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process  generating  the  water  surface  is  narrow-band  with  a  central  f re¬ 
al 

quencyjfl*  2.37  c/s.  The  water  waves  thus  generated  propagate  with  a 
velocity  (Stoker,  1957) 


v  ■  —  *  0.6584  m/s  (81) 

0) 

where  G  *  9.81  m/s2  is  the  gravitational  acceleration.  Beard's  measure¬ 
ments  (1967)  show  that  the  water  surface  can  be  considered  partially 
"frozen"  and  moving  along  the  water  tank  with  the  group  velocity,  which 
is  one-half  of  the  phase  velocity  (81)  for  deep  water  waves.  It  was 
also  established  that  the  surface  is  strongly  anisotropic;  i.e.,  the 
correlation  distance  in  the  direction  perpendicular  to  the  waterflow  (y) 
is  much  larger  than  in  the  direction  along  the  water  flow  (x),  so  that 
the  assumption  of  o’'  i-dimensional  roughness  is  to  some  extent  justified 
(assumption  D,  p.  ?,  assumes  an  infinite  correlation  distance  in  the  y 
direction) . 

The  correlation  function  envelope  corresponding  to  a  (t)  in  (75) 
was  found  to  be  very  close  to  an  exponential  curve  for  larger  T  [see  Fig¬ 
ure  8  in  (Beard,  1966)];  however,  for  t  0,  it  has  so  far  not  been  pos¬ 
sible  to  accurately  measure  B(t)  and  its  derivatives;  i.e.,  direct  mea¬ 
surement  of  the  mean  square  slope  and  the  correlation  distance  defined 
through  it  by  (9),  has  not  yet  been  performed.  On  the  other  hand,  indi¬ 
rect  inference  of  Tenv  from  the  measured  data,  i.e.,  by  (70),  is  not 
very  accurate,  since  B"(t)  is  not  known  for  x  -*■  0.  By  use  of  various 
methods  of  estimating  :env  from  the  measured  data,  it  was  found  that 
Tenv  is  of  the  order  of  1  meter. 

Beard's  experimental  results  (1965,  1967)  are  plotted  as  a  function 
of  the  variable 


o  siniii 
w 

=  g 


X 


(82) 
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The  variable  <  used  in  the  above  analysis  and  defined  by  (48)  is 
therefore  related  to  g  by 


/•<  3  - p -  g  (83) 

T  '2  COSlJ> 

env 


For  the  values  D  *  4.115  m  and  <p  =  22.9°  used  in  the  experiment, 
this  makes 


& 


9.89 

T 

env 


(84) 


As  tenv  is  of  the  order  of  1  m,  but  its  exact  value  *8  not  known, 
the  value  Tenv  *  0.989  m  was  assumed  so  that  (84)  reduces  to  the  numeri¬ 
cally  convenient  relation 


Sk  ■  10  g  (85) 

In  Beard's  experiment  (1965),  the  random  amplitudes  of  two  compo¬ 
nents  of  the  scattered  field  in  phase  quadrature  were  investigated.  The 
phase  of  these  components  with  respect  to  the  coherent  part  of  the  field 
could  be  varied  at  will  [which  was  not  the  case  in  the  experiment  by 
Clarke  and  Hendry  (1964)],  so  that  all  characteristics  of  the  variance 
ellipse  in  Figure  5,  namely  S0,  ax,  oy,  oel,  oe2  and  K2  =  <4ax/o2in 
could  be  directly  measured. 

The  data  measured  by  Beard  for  £0  as  a  function  of  ow  in  one  of  the 
experiments  (Beard,  1965,  Figure  29)  are  shown  in  Figure  8.  This  is 
compared  with  the  theoretically  derived  formula  (57)  where  k,  given  by 
(48),  is  converted  to  g  through  (85).  Equation  (57)  also  contains  the 
additive  constant  4>m  given  by  (33).  This  constant  cannot  be  accurately 
determined,  since  4>m  is  of  the  order  of  hundreds  of  wavelengths,  and  <j>m 
is  determined  by  its  fractional  part:  can  therefore  not  be  accurately 

determined  because,  among  other  reasons,  the  antennas  are  not  truly 
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point  sources  and  therefore  D  cannot  be  determined  with  the  accuracy  of 
a  small  fraction  of  a  wavelength.  However,  this  constant  is  additive 
and  does  not  change  with  the  water  roughness  ow,  so  that  it  simply  rep¬ 
resents  a  reference  phase  which  does  not  affect  the  function  £0(aw), 
except  in  shifting  its  plot  up  or  down  with  respect  to  the  origin.  In 
plotting  Figure  8,  the  additive  constant  4>m  was  chosen  so  as  to  give  the 
best  fit  with  the  experimental  data.  It  can  be  seen  that  the  theoreti¬ 
cal  curve  calculated  by  (57)  is  in  good  agreement  with  the  measured  val¬ 
ues  shown  in  Figure  8.  The  orientation  of  the  variance  ellipse,  as 
shown  by  both  the  theoretical  analysis  and  the  experimental  data,  ro¬ 
tates  with  varying  ow;  this  does  not  happen  for  plane-wave  illumination 
and  scatter;  i.e.,  when  transmitter  and  receiver  are  far  apart  and  high 
above  the  rough  surface  (Fraunhofer  scattering). 


Good  agreement  is  also  obtained  on  comparing  the  theoretical  formu¬ 
la  (62),  with  £  and  n  given  by  (49),  with  the  experimental  data  shown  in 
Figure  9.  The  experimental  data  are  those  reported  by  Beard  (1967,  Fig¬ 
ure  16).  The  ordinate  plots  o/<^,  where  is  the  amplitude  of  the  smooth 
surface  coherent  component.  The  latter  is  a  rather  complicated  function 
of  the  illuminated  area  and  the  gain  function  of  the  antennas  (Beckmann 
and  Spizzichino,  1963),  but  in  this  case  represents  simply  a  multiplica¬ 
tive  constant  which  can  be  absorbed  into  n  in  (62),  so  that  (62)  has 
been  plotted  in  Figure  9  with  the  multiplicative  constant  so  chosen  as 
to  give  the  best  fit  to  the  measured  data. 


The  extreme  values  o2^,  o^in  could  be  obtained  by  using  (37)  and 
substituting  (60),  (61)  for  0X2  and  0y2;  a  less  tedious  method  is  to 
calculate  K  first  and  then  to  use  the  equations 


o2  +  o2  ■  02, 

max  min 


o 

max 

_2 

°min 


K2 


(86) 


where  the  right  sides  are  known  from  (62)  and  (66)  respectively;  on 
solving  (86)  one  obtains 
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Fig.  9. 


RMS  scattered  amplitude.  Experimental  data  bv 
Beard  (1966),  theoretical  curve  computed  bv  (62) 
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.  10.  Extreme  values  of  the  variance  of  the  scattered 
field.  Experimental  data  by  Beard  (1965), 
theoretical  curve  computed  by  (87). 
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Fig.  11.  Ratio  of  extreme  variances  of  the  scattered  field. 
Experimental  data  by  Beard  (1966),  theoretical 
curve  computed  by  (66).  For  the  reason  of  the 
discrepancy,  see  last  paragraph  of  Section  4,  p.  21 
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and  these  values,  using  the  same  multiplicative  constant  as  in  Figure  9, 
are  plotted  against  the  experimental  data  (Beard,  1965,  Figure  34)  in 
Figure  10. 

Finally,  Figure  11  shows  a  comparison  of  K2  *  °max^°min  as  calcu“ 
lated  from  (66)  and  (68)  and  as  measured  by  Beard  (1967,  Figure  14). 
Although  for  large  (ow  sini|>)/X  both  theoretical  and  experimental  data 
approach  unity,  there  is  not  even  qualitative  agreement  between  theory 
and  experiment  for  ow  -*  0.  The  reason  for  this  discrepancy,  as  has  al¬ 
ready  been  mentioned  in  the  remarks  following  (74),  is  that  for  ow  -*■  0, 
the  variance  ellipse  shrinks  to  a  point;  the  ratio  K2  =  °max^°mJn  *s 
then  given  by  the  limit  of  an  indeterminate  expression  0/0,  which  can  be 
resolved  only  after  the  second  application  of  L'Hdpital's  rule,  yielding 
a  critically  sensitive  value,  which  can  obviously  not  be  reproduced  by  a 
theory  using  approximations  such  as  those  used  in  the  above  analysis 
(neglection  of  lateral  scattering,  asymptotic  methods  approximating  in¬ 
tegrals,  Newton's  method  approximating  the  location  of  the  saddle  point, 
etc. ) . 


7.  CONCLUSION 

The  theory  developed  in  Sections  2  through  5  leads  to  qualitative 
agreement  with,  and  provides  insight  into,  the  experimental  measurements 
of  the  phase-quadrature  components  scattered  by  a  random  rough  surface 
as  performed  by  Beard  (1965,  1967).  The  theory  takes  into  account  the 
finite  distance  between  transmitter  and  receiver;  i.e.,  the  sphericity 
of  the  incident  and  scattered  waves  by  using  a  hybrid  geometrical-physi¬ 
cal-optics  approach  and  leads  to  a  complicated  phase  distribution  of  the 
scattered  field.  The  latter  is  calculated  by  use  of  simplifying  assump¬ 
tions  and  approximate  integrations.  Unlike  the  case  of  plane-wave  illu¬ 
mination  and  reception  in  the  Fraunhofer  zone,  the  case  under  investigation 
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lends  to  asymmetrical  phase  distributions  even  if  the  rough  surface  is 
symmetrically  distributed,  and  to  a  rotation  of  the  variance  ellipse 
with  increasing  roughness.  Mathematically,  these  complications  are  re¬ 
flected  in  the  quantities  in  the  fact  that  the  ouantities  describing 
the  scattered  field  are  determined  not  by  the  exponential  functions 
appearing  in  plane-wave  scattering,  but  by  confluent  hypergeometric 
functions. 


32 


APPENDIX 

(l  3 

The  Confluent  Hypergeometric  Function  jF  I— ;  x 


The  quantities  of  interest  associated  with  the  scattered  field  are 
in  this  report  expressed  in  terms  of  the  confluent  hypergeometric  func¬ 
tion 

1  13  ill 

,  ,  2  2  *  2  x2  2  *  2  *  2  x3 

=  ±+  3  x-t-  3^2!  +  3  57  3| 

2  2  *  2  2  '  2  *  2 

-  1  +  —  x  +  —  x2  +  —  x3  +  x1*  + 

3  10  42  216  **• 

which  is  easily  programmed  on  a  digital  computer.  However,  for  quick 
information  computations,  a  short  table  of  y;  xj  is  included 

here,  since  tables  of  this  function  (Airey,  1926,  1927)  are  not  readily 
accessible. 
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TABLE  1-A 


X 

F  (i.  A-  x) 

1  1\ 2 *  2'  ) 

X 

F(l.  I. 

l  1 V  2 ’  2’ 

0.00 

1.00000 

1.1 

1.52757 

0.02 

1.00671 

1.2 

1.59700 

0.04 

1.01349 

1.3 

1.67129 

0.06 

1.02037 

1.4 

1.75083 

0.08 

1.02732 

1.5 

1.83603 

0.10 

1.03436 

1.6 

1.92736 

0.15 

1.05233 

1.7 

2.02531 

0.20 

1.07086 

1.8 

2.13041 

0.25 

1.08997 

1.9 

2.24325 

0.30 

1.10968 

2.0 

2.36445 

0.35 

1.13001 

2.2 

2.63478 

0.40 

1.15098 

2.4 

2.94764 

0.45 

1.17262 

2.6 

3.31041 

0.50 

1.19496 

2.8 

3.73183 

3.0 

4.22221 

0.55 

1.21801 

3.5 

5.83596 

0.60 

1.24818 

4.0 

8.22631 

0.65 

1.26638 

4.5 

11.7973 

0.70 

1.29175 

5.0 

17.1722 

0.75 

1.31796 

5.5 

25.3164 

0.80 

1.34503 

6.0 

37.7301 

0.85 

1.37300 

6.5 

56.7504 

0.90 

1.40191 

7.0 

86.0296 

0.95 

1.43178 

7.5 

131.289 

1.00 

1.46265 

8.0 

201.510 
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